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Abstract 

In this paper we discuss the quantum potential approach of Bohm 
in the context of quantum cosmological model. This approach makes 
it possible to convert the wavefunction of the universe to a set of 
equations describing the time evolution of the universe. Following 
Ashtekar et. al., we make use of quantum canonical transformation to 
cast a class of quantum cosmological models to a simple form in which 
they can be solved explicitly, and then we use the solutions do recover 
the time evolution. 
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1 Introduction 



The problem of the origin of the universe is perhaps as old as the mankind, 
but only in the recent years our understanding of the nature had become 
sufficiently deep to address it in the framework of physics. The theories 
we have in hands nowadays make it possible to trace the evolution of the 
universe starting from tiny fraction of a second after the mysterious Big Bang. 
This is, without doubts, one of the greatest achievements of the modern 
science. However, during these investigations it became more and more clear 
that our present knowledge is not sufficiently developed to answer questions 
concerning the Big Bang itself and what happened before. Actually this 
last question may happen to be not well posed, as there are indications that 
"before Big Bang" time did not exist at all. 

This circle of problems is the subject of a newly developed branch of 
physics which incorporates recent developments in high energy physics, quan- 
tum field theory, and quantum theory of gravity, and is called quantum cos- 
mology. Even if in the recent years there is growing activity in this field, 
there is a number of problems which not only are unsolved, but concerning 
which there is even no consensus, as to which of many possible questions are 
relevant. Let us list some of these problems. 

• It is natural to attack the problems of quantum cosmology by making 
use of some hybrid theory constructed from quantum field theory and 
the theory of gravity. Such a theory is certainly not known, but what- 
ever final shape it would take, it will provide us with the class of "wave- 
functions of the universe" as solutions. The problem is how should we 
interpret such wavefunctions. In quantum mechanics, a wavefunction 
provides us with the probability of finding a system in some configu- 
ration. But this is not really what we expect of quantum cosmology. 
The problem of major interest is how universe evolved and what is the 
class of initial conditions leading to the universe we observe. One can 
also pose a question as to whether there are any alternatives, i.e., are 
there any alternative consistent universes. 

• It is well known that the fundamental symmetry of general relativity is 
the general coordinate invariance. This leads, through the Dirac pro- 
cedure, to the quantum theory which has reparametrization invariance 
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built-in by default. This means, in particular, that the theory is not 
sensitive to any time reperametrizations, and therefore happens to be 
completely time-independent. Then how we can talk about 'evolution' ? 
This is the problem of time, which is probably the most important con- 
ceptual problem of quantum gravity |Q. 

• Equally well known is the fact that since the equations governing quan- 
tum mechanical evolution are linear, any linear combination of their 
solutions is a solution itself. But how could we interpret such a solu- 
tion consisting of, say, a sum of two different wave functions. Should 
we interpret this solution as describing two universes, or maybe hope 
that for some mysterious reasons, only one of them is relevant for our 
observations of the universe? This is the problem of decoherence (for 
review, see 0, section 5.5 and p| and references therein.) 

In the present paper we describe the quantum potential approach to quan- 
tum theories proposed by David Bohm || and developed in the context of 
quantum gravity in and for cosmology in |J and ||. This approach makes 
it possible not only to identify trajectories associated with 'wave function of 
the universe', but also to circumvent both the time and decoherence prob- 
lems mentioned above. By making use of the important technical result of 
Ashtekar, Tate, and Uggla 0, ||, who were able to transform the Hamilto- 
nian constraint of some interesting cosmological models to the form of free 
two or three dimensional wave equation, we are able to find a large class of 
exact wave functions for these models and to identify the corresponding time 
evolutions. 

Our investigations are undertaken within the framework of minisuper- 
space approximation to quantum gravity. This means that out of the infinite 
number of degrees of freedom of the full theory, we choose to investigate the 
dynamics of only finite number of modes. This approach is frequently crit- 
icized because there is no apparent reason why such a handicapped theory 
is to describe the evolution of real universe. However, one may argue that 
such a reduced model can result from taking into account the solutions of the 
full quantum theory with high degree of symmetry, exactly in the same way 
as one derives the Friedmann-Robertson-Walker cosmological model from 
classical general relativity. In fact, there is no apparent reason whatsoever 
for the FRW class of solutions of Einstein equations to describe the phys- 
ical universe, but surprisingly enough it does with great accuracy. There 
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is, no doubt, a great mystery in the fact that our (classical) universe is ho- 
mogeneous and isotropic, but there is no clear reason why this high degree 
of symmetry should not survive quantization. If the quantum universe also 
possesses a high degree of symmetry, the minisuperspace approach may be 
adequate after all. 

The plan of the paper is as follows. In section 2 we discuss the general 
features of the quantum potential approach. In section 3 we introduce the 
class of metrics which will be of our interest and derive the form of classical 
and quantum Hamiltonian constraints. In the next section, following ||, we 
analyze the canonical transformations which drastically simplify the form of 
Hamiltonian constraint, in fact in all considered cases this constraint reduces 
to the two- or three-dimensional wave equation. 

Having solutions of equations is, as it is well known, not sufficient. One 
must decide which of the solutions are 'physical'. The simplest criterion is 
to demand that the resulting wave function is normalizable. Therefore, in 
section 5 we address the question of what the inner product is. Section 6 
is devoted to description of a simple class of solutions of resulting equations 
and to the derivation of the time evolutions corresponding to the so obtained 
'wave functions of the universe.' 

2 The quantum potential 

In this section we will explain how the quantum potential approach works in 
the context of quantum cosmology. The general features of this approach (in 
particular, applied to the standard quantum mechanics and quantum field 
theory) are discussed in details in || and we recommend the reader to check 
these papers for general background and discussion of conceptual issues. 

For the cases which will be the object of our investigations in the sections 
to come, it is sufficient to consider a model for which the whole of quantum 
dynamics resides in the single equationF] 

H^z*) = Q^ViV,- - VV)) = Q° - VV)) M = 0, (1) 

what follows all quantum mechanical operators will be written in sans serif type 
face: a, b, . . . , A, B etc. 
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where g l i may be x- dependent. From now on we will call ip the wave function 
of the universe. Let us assume that ip has the following polar decomposition 

$ = Rix^exp (^(x 1 )) (2) 

with both R and S real. Inserting (0) into flU), we obtain two equations 
corresponding to real and imaginary part, respectively. These equations read 

Equation (^) will not concern us anymore. In this paper we assume that the 
wave function ip is a solution of equation (|l|), and thus this equation is iden- 
tically satisfied (even though we may not know what the explicit form of the 
wavefunction is.) On the other hand, equation ([3]) is of crucial importance. 
This equation can be used to derive the time dependence and then serves as 
the evolutionary equation in the formalism. 

For, let us introduce time t through the following equation 

t& „ 5H[S(x)\ 

dt 9 5{dS/dxi)' 1 ' 

This equation defines the trajectory x l (t) in terms of the phase of the wave- 
function S. Now we can substitute back equation (H) to @. Assuming that 
the matrix g l i has the inverse g^-, we find (x % 



•^i dx 1 

~ dt 



We see therefore that the quantum evolution differs from the classical one 
only by the presence of the quantum potential term 

h 2 l 

-V quant {x l ) = —-DR 
Z It 

on the right hand side of equation of motion. Since we assume that the wave 
function is known, the quantum potential term is known as well. 
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Equation (^|) is not in the form which is convenient for our further inves- 
tigations. To obtain the desired form we define classical momenta 



_ 5H[S(x)] _ , 
Pl ~ 5(dS/dx*) ~ 9ljX 

and cast equation (0) to the form 

H = \g ij PiPj + V{x*) - Y^ aR = °- ( ? ) 

We regard H, as the generator of dynamics acting through the Hamilton 
equations 

m 

% dx { 
dpi 

Thus we developed the scheme in which we can identify the time evolution 
corresponding to the wavefunction of the universe. This time evolution is 
governed by equations (§) subject to the constraint for initial conditions (|7|). 
This completes the technical part of quantum potential program. However a 
number of remarks is in order. 

1. Equation (|8|) with the hamiltonian TC defined by (|7|) is identical with 
the classical equations of motion with the only difference that in the 
hamiltonian the classical potential is appended by the quantum term 
y quant which is a quantity of order h 2 . Observe that in the course of 
deriving the quantum potential equations of motion, we did not make 
any approximations and therefore the resulting theory is completely 
equivalent to the original quantum system. What makes this approach 
different from the standard quantum mechanics, is that the wave func- 
tion whose interpretation is rather obscured, especially in the case of 
quantum gravity and quantum cosmology, is now replaced by trajecto- 
ries with well understood physical meaning. 

2. The quantum potential interpretation may be used to obtain a well 
defined semi-classical approximation to quantum theory. Indeed, it 
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can be said that the system enters the (semi-) classical regime if the 
quantum potential is much smaller than the regular potential term. 
This observation has been extensively used in the paper [||, where we 
were interested in the various mixed regimes (quantum matter and 
classical gravity etc.) 

3. One of the major adventages of the quantum potential approach is that 
it provides one with an affective and simple way of introducing time 
even if the system under consideration has a hamiltonian as one of the 
constraints. In particular this approach serves as a possible route to 
final understanding of the problem of time in quantum gravity. We do 
not intend to suggest that this is the only possible approach to this 
problem, however, to our taste, this is the best one available now. 

4. Related to this is the problem of interpretation of wave functions which 
are real. This problem has been a subject of numerous investigations, 
but from the point of view of quantum potential the resolution of it 
is quite simple. We just say that real wavefunctions (of the universe) 
represent a model without time evolution (and therefore time) at all. 
Recall that the standard (time-dependent) Schrodinger equation does 
not have any real solutions. It is clear from the formalism: x is just 
equal to zero, so nothing evolves and therefore there is no clock to 
measure time. On the other hand, equation (^) means that for the real 
wave function the system settles down to the configuration for which 
the total potential (i.e., classical plus quantum) is equal to zero. 

5. There is one particular example of equation (1) which will be important 
for us later on. Consider the situation when the classical potential is 
absent. Since the kinetic term □ is a hyperbolic differential operator, 
the wave function is a sum of function corresponding to 'retarded' and 
'advanced' waves, to wit 

■ip = ipi(u) + ip(v). 

Now the quantum potential is defined in terms of the same operator 
□ acting on the modulus of the total wave function. We see therefore 
that if we want the quantum potential not to vanish, we must keep 
both components of the wavefunction decomposition above. 
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6. The definition of time by equation (^) is, of course, not unique. In fact, 
we can can use a more general expression 

V = N(x) 6H[S{X)] (9) 

where, in the case of gravity, the function N is to be identified with 
the lapse function of ADM formalism. 

This completes our formal investigations of the quantum potential ap- 
proach. Let us now turn to cosmology. 



3 Classical cosmological models 

This section concerns a class of homogeneous cosmological models, called the 
diagonal, intristically multiply transitive models, whose quantum evolution 
we will investigate in the following sections. Our presentation essentially 
follows the paper || and we recommend the reader to consult this paper for 
more details and references. 

A spacially homogeneous 4-metric can be expressed in the form 

ds 2 = -N\t)dt 2 + j2gu{u l )\ (10) 
i=i 

where N(t) is the lapse function and uo l are some appropriate 1-forms on the 
3-space. To make this metric compatible with the non-evolutionary Einstein 
equations, to 1 must satisfy certain conditions whose explicit form will not 
concern us. Particular forms of id 1 are classified and called Bianchi types. 

Since the metric coefficients gu are positive, they can be represented in 
the form 

9u = e p 

and it turns out that the following parametrization, due to Misner, is very 
convenient 

/ P 1 \ /ll \/3 \ / 0° \ 

P 2 = 1 1 -V3 \\ P + \, 
\(3 3 J V 1 "2 )\(3- ) 
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(x°, x + , x-) = -^(2/3° - -p° + 2/3+, V3/T). 

Defining the momenta associated with the variables x l 

{x\ Pj } = 5} 

and turning to the so called Taub time gauge which fixes the lapse function 
N to be 

N T = 12exp(3/3°) 

we finally arrive at the following simple form of the single remaining Einstein 
equation (in phase space) 

H = H + H + + H- = 0, (11) 
Ho = -\ P l + k e^°, (12) 

H + = ipl + k + e-^ + , (13) 

H- = \k_pl (14) 

where the coefficients ko, k± for various models Bianchi are given in the table 
above. 

Before going any further let us make an important remark. As it is well 
known the Einstein lagrangian for cosmological models has the generic form 



^ J dt(px-N(t)H). 



Setting M P = 1, as it is usually done, means that all quantities are scaled 
by the Planck mass, length etc. Thus if in the numerical computations some 
quantity is equal 1, this quantity is 1 times an appropriate power of M P 
and this means that 1 is just the border between the "classical" and the 
"quantum" world. 

There is one more model, the so called Bianchi type V model for which 
the scheme is applicable as well; for this model, the x and (5 variables are 
identical and the coefficients are given by k = 72, k + = 0, k^ = 1. 

All the models listed above have the important property of being separa- 
ble which makes it possible to find a canonical transformation leading to the 
potential-free theory. 
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Table 1: Coefficients of hamiltonian constraint for various Bianchi models 

4 Classical and quantum canonical transfor- 
mations 

It is well known that canonical transformations are a very convenient tool 
of classical mechanics. If applicable, they are used to simplify complicated 
mechanical systems to the ones which are easy to deal with. It turns out 
however that in the case of quantum mechanics our understanding of canon- 
ical transformations is surprisingly poor. To our knowledge the only strong 
result is due to Berezin |§ and covers only linear canonical transformations. 
Therefore, in a non-linear case we cannot refer to any theory, but rather we 
must just check if a classical canonical transformation can be elevated to the 
quantum regime. The problem is that we must make sure that all the pe- 
culiarities of quantum theory (like, for example, operator ordering problem) 
do not spoil the classical transformation. This is exactly what we will do in 
this section in the context of models described in section 3. 

The generic piece of the hamiltonian constraint can be written as 

e-(p 2 + aexp bx), 

where e = ±1, and a, b are coefficients which can be read off from the table 
in the previous section. The obvious goal will be therefore to replace this 
term, after canonical trasformation, by a square of pure momentum. Thus 
the canonical transformation must read 

P = ^/p 2 + aexp 2bx. (15) 

To define the new variable X we must solve the equation {Q, P} = 1, Now 
we must distinguish two cases depending on the sign of the coefficicent a 
above. 

X = i (log[-p + vV + ae 2bx ] - logtvW"]) , a > 0, (16) 
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X = X - (\og\ V - sjp 2 + ae 2bx ] - \og[V^e bx ]^j , a < 0. (17) 

The inverse canonical transformation reads 

p 

^ae bx = — — — , p= -Ptanh&X, a > (18) 
cosh bX 

V^ae bx = — — — , p = PcothbX, a < 0. (19) 
sinh bX 

Let us observe that the momentum P defined by (|T5|) is positive. This does 
not pose any problems in classical theory, but we will have to struggle a 
little with this condition in the quantum case below. In the paper || one 
can find a detailed discussion of topology of the phase space of the models 
after canonical transformation. It should be stressed that since the canonical 
transformations lead to a free theory, this topology encodes all the nontrivial 
information concerning the dynamics of the model. 

Let us now turn to the quantum picture. As we mentioned above there is 
no complete theory of quantum canonical transformations, and all particular 
(nonlinear) cases must be discussed separately. In what follows, we will 
discuss the case of positive a, the negative a case can be discussed similarly. 

We look for the linear transformation which maps wave functions ip{x) 
into wave functions ty(X). In fact the canonical transformation in our 
case is nothing but the transformation to the basis in which the operators 
p? + <2j exp bi Xj are diagonal. It follows from the basic principles of quan- 
tum mechanics that such a (unitary and linear) transformation exists. This 
transformation must have the property that solutions of constraint equations 
before and after the canonical transformation are mapped into each other. 
Moreover we would like this transformation to be unitary; this, however re- 
quires the knowledge of the inner products. Since the inner product in the 
free case is easy to construct we use the transformation to define the inner 
product in the space of original wave functions. 

Thus we start with writing 

/+oo 
dX < x\X > *(X). 
-oo 

To find the kernel < x\X > satisfying the requirement above we investigate 



the action of the quantum operators corresponding to flT8|) on both sides of 
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the above equation. To this end, we must choose some particular operator 
ordering, and it turns out that only the ordering "PQ" leads to the desired 
result. 

We therefore assume that 

/+oo 
dX < x\X > PtanhoX^(X), (20) 
-oo 

/+oo 1 
dX < x\X > P — ——*(X), (21) 
-oo COSh OA 

From these two conditions we have (in position representation) after inte- 
grating by parts 

d d 
-i— < x\X >= -2ttt- (< x\X >) tanh&X, 
ox oX 

and 

d 1 

y/aexp(bx) < x\X >= i-^rp (< x\X >) 



'dX K 1 cosh&X' 
which can be solved to give 

< x\X >= exp ^-^e 6x sinh(6X)^J (22) 

In general two or three dimensional case the kernel has the form 

/ 3 



< xi,x 2 , x 3 \Xx, X 2 , X 3 >= exp [ -i -^-e^ 1 sinh(6JQ 



i=i 



It should be observed that in order that the procedure described above make 
sense and that solutions of constraints are mapped to solutions, the wave 
function ^(X) must be integrable with the kernel fl22"|) at least together with 
its first derivative. By changing variables in the integral sinh bX = Y, we see 
that 

*(,) = £V exp {-if e-r) 1 ^ f ,nX(Y)) 

and the integral exists (in the sense of principal value) if ^(Y) tends to a con- 
stant when Y — > ±oo. Of course, ^(Y) cannot have any nasty singularities 
in finite interval. 
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Therefore given any solution of simple free system we can, in principle, 
construct a solution of the original theory. The problem is that the integral 
transform defined above is very complex because the kernel oscillates very 
rapidly for large X. For that reason it is not only hopeless to find the 
integral in terms of tabulated functions but it is even difficult to compute 
the transform numerically. However some approximate methods (for large x, 
for example) may work for some ^(X). 

5 The inner product 

We succeeded therefore to reduce our problem to the problem of free scalar 
particle in two or three dimensions. However not all solutions of the con- 
straint equation 

nV(X) = (23) 

can be interpreted as a physical wavefunction. First of all, in order to make 
contact with the original theory, we must make sure that the integral trans- 
form defined in the previous section exists. Secondly, the wavefunction is 
supposed to describe physical probabilities, thus it must be normalizable. 
Thus we must define an inner product in the space of solutions of equation 

The problem of defining the inner product for Dirac quantization scheme 
is unsolved in general (see however fTOfl). In the case in hand, however, it 
suffices just to construct the inner product. 

To this end, let us observe that any solution of equation 
momentum representation, written in the form 

*(P) = StftPiPjfiiP), 

where ^f(P) is an arbitrary function of P sufficiently well behaving at rj^PiPj = 
0. In any other representation one must replace ^{rfi PiPj) with the delta 
function of the corresponding operator 

/+oo 
dX exp(iAj7 y PiPj-)tf. (24) 
-oo 

Now we can define the inner product on the space of solutions of costraint 
equation to be 

|| || 2 = [ ^StfiPiPj)*. 



23) can be, in 
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This inner product is not very convenient for our future investigations. The 
reason is that the wave function being a solution of (two dimensional) wave 
operator has a general form 



/+oo r+oo 
dkA k e lkU + / d\B^ w 
-oo J — oo 



where U = Xq + X% and V = Xq — X\ . For that reason we would like to have 
an equivalent definition of the inner product for wave functions in position 
representation. 

Following |Hj we proceed as follows. We define the inner product to be 



-oo 

j2 



II y{x) || = / d 2 x^*(x)fi^(x), 

J — oo 

where fi is a gauge fixing operator. In the case in hand this operator is chosen 
to be |p^-Q This operator acts in a simple way in momentum representation. 
We have therefore]] 

ii ir= 

= / d 2 X'd 2 Pd\V*(X') < X'\P > exp -A — < P\X > V(X) -- 

J-oo \ 2P oPq J 

= / +00 d 2 X'd 2 PdA^(XOe- i ( x ^o+xiP 1 ) e i((Xo(Po-^)+x 1 p 1 )^ (Xo)Xi) 

J — oo 

Integrating over P% and X[, we get 

dX' Q dX x dX dP Q d\^{X^ X 1 )e- iX '° Po e i{Xo{Po ^ ) ^{X , X x ) = 



^ dX' Q dX 1 dX dP y*(X' ,X 1 )e iP °( x °- x o)5 (^J *(X ,X 1 
which after integration over Xq gives 

30 

dX' dXidP V*{X' , X 1 )e- lPoX o2P ^(0, X x ) = 



2 In general case when the comutator of gauge fixing operator and gauge constraint is 
not constant, the formula above must be generalized, see |1C{ |, 

3 < | > is the inner product of the original Hilbert (or, more generally Krein) space of 
the model. Thus < P\X >= e iPX . 
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/ + 00 (J 

This inner product can be also written as 



x' =o 



+00 



dXdY V*(Y,X) d Y V(Y,X) 



Y=0 



(25) 



This is in fact the standard inner product for massless spin zero particle, as 
it could have been expected. Observe that we clearly have here normalizable 
negative norm states. This fact is well known, and there is an open problem 
what is the physical meaning of such states, and/or if their existance indicate 
the need of third quantization. 



6 A simple model 

In this section we make use of the machinery built above to discuss a simplest 
possible nontrivial model, namely 

tt(X , X x ) = e i{k+l)u + e i{k - l)v , (26) 

where U, V — Xq ± X + as before. The above wave function is certainly not 
physical (it is not normalizable), however it has the virtue to analogous to 
the plane wave states in the standard quantum field theory. For that reason 
all the results below should be taken with the grain of salt, however it does 
not seem unlikely that some physical prediction of this simple model may 
hold for the full theory which will be discussed elsewhere. 

Given the wavefunction ( p6|) , one can readily find the quantum potential. 
Equation (7) reads in our case 

-P 2 + Pi + (k 2 - I 2 ) = 

and from the hamiltonian equation of motion we find easily that 

P = k X° = -kt (27) 
P + = l X + = lt + t (28) 

From the condition for the canonical transformation to be meaningful, we 
see that both k and I must be positive. Now we are ready to turn to the 
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starting variables. To this end we choose to work with the Bianchi type IX 



model, and using equation (jTSf) , we find 



VA8e 2 ^ x ° = ^-7^, Po = -/tanh2v / 3H, (29) 

cosh 2v 3kt 

^~*" = cosh4^ + f ,) - P + = »^4V3( ii + t0 ). (30) 

Looking at the Misner parametrization, we see that the metric depends on t 
through the combination 

o _ 1 / k /2yg\ 1/4 co^^/m + h ) 
3 ° g ^4v^ V 1 J cosh(2 v / 3H) 

„. 1 , ( k cosh(4\/3^ + t ) \ 

P = r ^g 



6 *\2l cosh(2V3A;t) / 

Now we can discuss the characteristic features of the model. It is described 
by the metric 

ds 2 = -N 2 dt 2 + e 2/3 ° (e 2/3 V) 2 + e 2/3 V) 2 + e~ 4/3 V) 2 ) , 

where e 2f3 ° is the scale factor and e 2l3+ is a measure of the anisotropy. The 
physical time is given by the monotoneous function of the parameter t, 

cosh 1/4 (4 v / 3/t + t ) 



cosh(2 v / 3H) 



dt 



and for large t, r behaves as r ~ j^k^ 1 ^ ~ 2k \ The scale factor is a sym- 
metric function of t with one extremum, which is a minimum for I > 2k, 
and a maximum for I < 2k. Similarly, the anisotropy factor has only one 
maximum for 21 > k and only one minimum for 21 < k. Thus we have four 
qualitatively different classes of models, depending on different values of k, I. 
Observe that the classical behaviour corresponds to k = I (in this case the 
quantum potential vanishes.) 

The model above should be viewed only as an illustration of application 
of our method to quantum cosmological models. As it was mentioned at the 
beginning of this section its physical value is diminished by the fact that the 
corresponding wavefunction is not normalizable. The more realistic models 
will be considered in the separate paper. 
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